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Abstract. By reversing the time variable we derive a stochastic representation for back- 
ward incompressible Navier-Stokes equations in terms of stochastic Lagrangian paths, 
which is similar to Constantin and Iyer's forward formulations in [5] . Using this represen- 
tation, a self-contained proof of local existence of solutions in Sobolev spaces are provided 
for incompressible Navier-Stokes equations in the whole space. In two dimensions or large 
\ viscosity, an alternative proof to the global existence is also given. Moreover, a large de- 

Ph ■ viation estimate for stochastic particle trajectories is presented when the viscosity tends 

^ (— I to zero, 

-i— > 

B: 

1. Introduction 

The classical Navier-Stokes equations describe the evolution of velocity fields of an 
incompressible fluid, and takes the following form with the external force zero: 

[ d t u + (u ■ V)u - Mm + Vp = 0, t ^ 0, 

d 



00 

o 



V-tt = 0, u{0) = u , ^ 

where column vector field u = (u 1 , u 2 , u 3 )* denotes the velocity field, p is the pressure and 
v is the kinematic viscosity. When the viscosity v vanishes, the above equation becomes 
the classical Euler equation: 



fd t u + {u- V)u + Vp = 0, 0, 

\V-m = 0, u(0) = u 0} ^ ' 

which describes the motion of an ideal incompressible fluid. The mathematical theory 
about Navier-Stokes equations and Euler equations has been extensively studied and the 
existence of regularity solutions is still a big open problem in modern PDEs. 

Recently, Constantin and Iyer [6] presented an elegant stochastic representation for in- 
compressible Navier-Stokes equations based on stochastic particle paths, which is realized 
by an implicit stochastic differential equation: the drift term is computed as the expected 
value of an expression involving the stochastic flow defined by itself. More precisely, let 
(u, X) solve the following stochastic system: 

X t (x) =x+ f u s {X s {x))ds + V2^B t , t ^ 0, 
[u t = EP[(V t A t )(u oA t )], 



Keywords: Backward Navier-Stokes Equation, Stochastic Representation, Global Existence, Large 
Deviation. 
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where B t is a 3-dimensional Brownian motion, P is the Leray-Hodge projection onto 
divergence free vector fields, and VM t denotes the transpose of Jacobi matrix VA t . Then 
u satisfies equation ([T]) with initial data Uq. One of the proofs given in [B] is based 
on a stochastic partial differential equation satisfied by the inverse flow A. By using this 
representation, a self-contained proof of the existence of local smooth solutions is provided 
in [HI [13]. 

Let (u,p) solve ([1]). Notice that if we make the time change: 

u(t, x) := — u{— t, x), p(t, x) = p(—t, x) for t ^ 0, 
then u satisfies the following equation (called backward Navier-Stokes equation here): 

r d t u + (u ■ V)m + vAu + Vp = 0, t ^ 0, 

\V-fi = 0, u(0) = u . ^ 

The purpose of the present paper is to give a slightly different representation for u by 
using backward particle paths. More precisely, let (u, X) solve the following stochastic 
system 

X t>s (x) = x + [ u r (X t , r (x))dr + V2^(W S - W t ), t ^ s ^ 0, 

Jt (5) 

u t = ^P[(y t X t>0 )(u oX tfi )], t^O. 

Then u satisfies backward incompressible Navier-Stokes equation (jl]) with final value Uq. 
Intuitively, by reversing the time, the starting point is changed as the end point. Hence, 
representation ([5]) is essentially the same as Q. But, equation ([5]) is easier to be dealt 
with mathematically since the unpleasant term A = X -1 which usually incurs extra 
mathematical calculations does not appear in Such a representation will be proved 
in Section 2. We emphasize that representations ([3]) and §5§ are useful in numerical 
computations (cf. [EKE])- Direct calculations shows that the second equation in §5§ is 
equivalent to 

Qt = V x ut = E[(V- 1 X tj0 )(o;o o X tj0 )}, (6) 
ut = -A- : V x u t , (7) 

where uj t is the vorticity, V _1 X ti o is the inverse of Jacobi matrix VX^o and (J7|) is exactly 
the Biot-Savart law. 

We also mention other stochastic formulations for incompressible Navier-Stokes equa- 
tions. In [5], a representation formula for the vorticity of three dimensional Navier- 
Stokes equations was given by using stochastic Largrangian paths, however, there is no 
a self-contained proof of the existence given there. In [TS], Le Jan and Sznitman used a 
backward-in-time branching process in Fourier space to express the velocity field of a three- 
dimensional viscous fluid as the average of a stochastic process, which then leads to a new 
existence theorem. In [3], basing on Girsanov's transformation and Bismut-Elworty-Li's 
formula, Busnello introduced a purely probabilistic treatment to the existence of a unique 
global solution for two dimensional Navier-Stokes equations, where the stretching term 
disappears, and the non-linear equation obeyed by the vorticity has the form of Fokker- 
Planck equation. Later on, Busnello, Flandoli and Romito in [1] carefully analyzed an 
implicit probabilistic representation for the vorticity of three dimensional Navier-Stokes 
equations, and a local existence was given. In that paper, much attentions were also paid 
on a probabilistic representation formula for a general system of linear parabolic equa- 
tions. Moreover, in [31 H], an interesting probabilistic representation for the Biot-Savart 
law was also given and analyzed so that they can recover the velocity from the vorticity 
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by probabilistic approach. Recently, Cipriano and Cruzeiro in [5] described a stochastic 
variational principle for two dimensional incompressible Navier-Stokes equations by using 
the Brownian motions on the group of homeomorphisms on the torus. More recently, 
Cruzeiro and Shamarova [7] established a connection between equation (TjJ and a system 
of infinite dimensional forward-backward stochastic differential equations on the group of 
volume-preserving diffeomorphisms of a flat torus. 

This paper is organized as follows: In Section 3, we shall give a self-contained proof 
of local existence in Sobolev spaces. The proof is based on successive approximation 
or fixed point method as in [12]. Therein, Iyer considered the spatially periodic case 
and worked in Holder continuous function spaces. When Sobolev spaces are considered, 
we have to overcome the difficulty due to the non-closedness of Sobolev spaces under 
pointwise multiplications and compositions. Thus, it seems to be hard to exhibit the 
same proof in Section 3 for representation ([3]) due to the presence of A = X" 1 . A key 
point in the proof lies in that the flow map x \— > X t s (x) preserves the Lebesgue measure, 
i.e., det(VX M ) = 1. 

In Section 4, we shall give an alternative proof to the global existence when the spatial 
dimension is two or the viscosity is large enough in any dimensions. Such results are well 
known. For two dimensional Navier-Stokes and Euler equations in the whole space, the 
global existence of smooth solutions are referred to [19] . The global existence of regularity 
solutions for large viscosity is referred to [IB]. Recently, Iyer [T4] presented an alternative 
proof to the global existence for small Reynolds. His proof is based on the decay of heat 
flows and stochastic representation ([3]). Following [14], we will give a different proof based 
on Bismut's formula and representation ([S]). 

Let (u u ,X u ) denote the solution of equation (jSJ). In Section 5, as v goes to zero, an 
asymptotic probability estimate of X u in diffeomorphism group is presented by the well 
known large deviation estimate for stochastic diffeomorphism flows. 



2. Stochastic Representation of Backward Incompressible Navier-Stokes 

Equations 

We begin with some notational conventions. Fix d ^ 2 and put 

N := {0} U N, R_ := (-oo, 0], I := (d x d)-unit matrix 

and 



I 2 



Mlja 



< +oo 



For a differentiable transformation X of M. d , the Jacobi matrix of X is given by 



VX 



/ <9iX\ d 2 X\ 

d 1 x 2 , d 2 x\ 

\ diX d } d 2 X d , 



ddX 1 \ 
d d X* 

d d X d j 



where d t = We use V*X to denote the transpose of VX. For k <G N , let Cj;(R d ] R d ] 
denote the space of fc-order continuous differentiable vector fields on M d with the norm: 



u 



sup \D a u{x)\ < +oo, 



where D a denotes the derivative with respect to the multi index a. 
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Let u G C(R_; Cf R d )) and t ^ a t = a(t) G Lf oc (R^; R d x Z 2 ) satisfy 

oo 

J> fc (tK(t)=I. (8) 

fc=i 

Let {Xt, a (x),t ^ s < 0} solve the following SDE 

PS PS 

X t!S (x) = x + u r (X t ,r(x))dr + V2^ J a r dB r , (9) 

where B t := {B>°, KO.fcGN} is a sequence of independent standard Brownian motions 
on some probability space P). Thanks to ([8]), the diffusion operator associated to 

equation ([9]) is given by 

L t g = uAg + (u t ■ V)g. 
We first prove the following result. 



Theorem 2.1. Let 



be a C 2 -vector field satisfying 



\D a <f>(x)\^C{l + \xf), |«K2, (3>0, 
and f G C(1R_; C%(R d ; R d )). Define 

A^(t,x) := (V% o (x))0(X t)O (x)), 



A/(*,aO 



(V%, r (x))/ r (X 4 , r (o;))dr 



and 



:= EA^t, x) - EA f (t, x). 
Then w G C 1,2 ((— oo,0) x R ') satisfies the following backward Kolmogorov's equation: 

d t w t + L t w t + (V^tjwt = ft, limw t (x) = <f)(x). (10) 

Proof. Let g be a twice continuously differentiable function satisfying 

|£> Q £(x)| < C(l + Ixl' 3 ), |a|<2, /3>0. 
For /i > 0, by Ito's formula we have 



Eg{X t ^ t (x)) =g(x)+E 



(L r g)(X t - hjr (x))dr 



t-h 



From this, it is easy to see that 
1 



hi 



Eg(X t - h>t (x)) - g(x) 



(L r g')(Xt_/ l)r (a;))dr 



— > (L t g)(x) as /i — > 0. 



(11) 



Noticing that 
we have 

VI t _ M (i) = (V^ )o^ M (i) • VI(_ M (x) 
Thus, by the independence of Xt-h,t(x) with X ti0 (x), we have 

EA^t -h,x) = E[(V%_ M (a;))A^t, X t _ M (z))] 



E 



(V t X t _ M (x))((EA^(t)) o X t - htt (x)) 



and 



EA 7 (t-M) = E [(V t X t _ M (x))A / (t, X t _ M (a;))] 

4 



+E 



(V%_ v (x))/(r,X 4 _,, r (x))dr 



t-h 



E [(V%_ M (x))((EA / (i))oI t _ M (i))] 



-E 



(V t X t -. htr (x))f r (X t - h , r (x))dr 



t-h 



Hence, we may write 

w t -h{x) 



E 



{V t X t . htt {x))w t (X t . ht t(x)) 



-E 



(V%_ v (x))/ r (X t „ v (x))dr 



t-h 



and 



- w t -h{x)) 



4- 

l r 

~h 

+ yE 

n 

rh 



V'Xt^tix) -I)w t (X t ^ h , t (x)) 
E (w t {X t - ht t(x))) - w t (x) 



(V t X t _ /l , r (a;))/ r (X t „ v (x))dr 



t-h 



lZ(t,x) + I»(t,x)+lS(t,x). 



Observing that 



we deduce 



VX t u t (x)-l 



(Vw s ) o X t - hjS (x) ■ VX t _ hjS (x)ds. 



t-h 



By (llip we have 



\jmI*(t,x) = -[(V t u t )w t ){x) 



liml%(t,x) = -(L t w t )(x). 

h|0 



Moreover, a simple limit procedure also gives 

\wal%(t,x) = f t (x). 

Combining the above calculations, we conclude that 



]im~(w t (x) - w t -h{x)) = -(L t w t )(x) - [(V'uj)^^) + f t (x) 

hlO tl 

Equation (ITU]) now follows (see (TTJ p. 124] for more details). 



□ 



Remark 2.2. A more general Feynman-Kac formula for a deterministic system of para- 
bolic equations was given in [I] . However, the proof is simpler in our case. In representa- 
tion (T3j), if we define w t := E[(VMf)(w ° A t )], then w t also satisfies ( Q7| ) with f = 0(see 
@ p.343, (4.5)];. 

Basing on this theorem, we can give a stochastic representation for backward Navier- 
Stokes equation (jl]) as in [6]. 
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Theorem 2.3. Let v ^ anduo G C^(M. d ; M. d ) a deterministic divergence- free vector field, 
and f G C(R_; C%(R. d ; M, d )). Suppose that a satisfies (TJJ), and (u,X) solves the stochastic 
system: 

X tiS (x) = x + I u r (X t , r (x))dr + V2u f a r dB r , t ^ s ^ 0, (12) 
Jt Jt 

u t = PEA^(t)-PEA^), t^O, (13) 

where P is the Leray-Hodge projection onto divergence free vector fields, and A^ o and 
are given by 

AMt.z) := (V%,o(x))«o(X t , (x)), 



A?(*,z) := y (V% >r (x))/ r pf tjr (x))dr. 

T/ien it satisfies the backward incompressible Navier-Stokes equation: 

d t u + (u ■ V)u + z/A?i + Vp = f, t ^ 0, 
V • it = 0, w(0, x) = Uq{x). ^ 
Conversely, if u solves backward Navier-Stokes equation [L~f\ ), then u is given by (flff|). 
Proof. First of all, let 

w t (x):=EAl o (t,x)-EA u f (t,x). (15) 
By Theorem 12.11 w(t,x) = w t (x) solves the following backward Kolmogorov's equation: 

d t w + (u ■ V)w + (V%)w + uAw = /, tu(0, x) = u (x). (16) 
In view of w = Pw, we may write 

w = u + Vg. 
Substituting it into equation (fTBl) . one finds that 

c^w + (it ■ V)m + j^Am + Vp = /, m(0, x) = uo(x), 
where ^ 

p = d t q + (w • V)g + z/Ag + — | | 2 . 

Conversely, let (u,p) solve ([II]) . As above, if we define u> by (ITBI) , then to satisfies 
equation (ITBI) . We need to show that it = Pw, or equivalently, for some scalar valued 
function q 

v := w — u = Vg. 
By (fTBj) and (fbil) . t> solves the following equation: 

d t v + (u ■ V)v + (V*u)i; + z^Ai; = Vp - -V|u| 2 , w(0,x) = 0. (17) 



Let 

g(t,x) :=E^jT [^\u(r,X t;r (x))\ 2 - p(r,X t>r (x))]dr 
Then g solves the following equation (cf. |llj): 

d t q + (u ■ V)g + uAq = p — - |w| 2 , g(0, x) = 0. 
Taking gradients for both sides of the above equation yields 

d t Vq + (u ■ V)Vg + (V*u)(Vg) + z/AVg = Vp - ^V|w| 2 . 
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By the uniqueness of solutions to linear equation ffTTj) . v = Vg. □ 

Remark 2.4. In the above proof, we have assumed that the solutions are regular enough 
so that all the calculations are valid. The existence of regular solutions will be proven in 
the next section. 

3. A Proof of Local Existence in Sobolev Spaces 

With a little abuse of notations, in this and next sections we use p to denote the 
integrability index since the pressure will not appear below. For k £ No and p > 1, 
let W k ' p {R d ;R d ) be the usual Revalued Sobolev space on R d , i.e, the completion of 
C£°(R d ; R d ) with respect to the norm: 

k 

\\ u \\k, P ■= \\u\\ p + 22 II VJm IIp> 

where || • || p is the usual Z^-norm, and V- 7 is the j-order gradient operator. Note that 
W°> p (R d ; R d ) = L p (R d ; R d ) and the following Sobolev's embedding holds: for p > d (cf. 

HDD 

W 1,p (R d ;R d ) ^ L°°(R d ; R d ), i.e. \\ ■ (U ^ c|| • \\ liP , (18) 

where c = c(p, d) . Below, we shall use c to denote a constant which may change in different 
occasions, and whose dependence on parameters can be traced carefully from the context. 
Let W^ p (R d ; R d ) be the local Sobolev space on R . We introduce the following Banach 
space of transformations of R d : 

X k+2, P ._ f x e w k+2,p( Rd . Rd ) . + nvxiU + ||V 2 X|| fciP < +oo|. 

Definition 3.1. The Weber operator W : L p (R d ; R d ) x X 2 ' p -> L p (R d ; R d ) is defined by 

W(vJ) :=P[(V^)v], 

where P is the Leray-Hodge projection onto divergence free vector fields. 

Remark 3.2. P = I — V(— A) _1 div is a singular integral operator(SIO) which is bounded 
in L p -space forp £ (1, oo) (cf. [22] ). 

We now prepare several lemmas for later use. 

Lemma 3.3. (i) For any k £ No and p > d, there exists a constant c = c(k,p,d) > 
such that for all v £ W k+2 ' p (R d ; R d ) and £ £ X k+2 ' p , 

\\VW(v,£)\\ k+1>p ^ cdlV^IU + \\V 2 i\\ ktP ) ■ \\Vv\\ k+ljP (19) 

(ii) For p > d, there exists a constant c = c(p,d) > such that for all Vi,v 2 £ 
W 2,p (R d ; R d ) and £ u £ 2 £ X 2 > p with l x - t 2 £ L p (R d ; R d ) 

\\W(v h £ 1 )-W(v 2 ,£ 2 )\\ p ^ c(\\v 1 \\ 2 J£ 1 -£ 2 \\ p +\\V£ 2 \\oo\\v 1 -v 2 \\ p y (20) 

Proof, (i) Noting that 

P((V*£)v) + P((VV)£) = P(V(£ ■ v)) = 0, (21) 

we have 

d t W(v,£) = P[(V t 9 i £)t» + {VH)div\ = V[-(V t v)d l £ + (W)<9^] 

and 

djdiW(v,e) = P[-(V t v)d j d i i - (V'd^dd + (VV)^-^ + (V^ftu]. 
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Hence, by f fTBl) we have 

||<9iWM)|| P «S c(||(VM^II P + ll(VV)^|| p ) 

< c(||V«|| p -||^||oo + ||V^|oo||9i«y 
*c c\\Vv\\ p - IIWHoo 

and 

WdjdiWivJ)^ ^ c(||V 2 £|| p - II VwHoo + HV^Hoo - II V 2 w||p) 
^ c(|| V 2 £\\ p ■ || Vv\\ 1>p + llV^Hoo - ||V 2 v|| p ) 
^ c(||W|U + ||V 2 £|| P )- || Vv\\ liP , 

which produces 

||vwm)|| 1jP ^ c(||v£|U + ||v 2 £|| p ) • \\v v \\ hp . 

The higher derivatives can be estimated similarly, 
(ii) By (JZED, we have 

= -P((VS)(*i - £ 2 )) + P((V t 4)(^i - v 2 )). 

So, 

|| W( Vl , h) - W(v 2 , l 2 ) \\ p ^ c\\ (V%) (£ 1 - l 2 ) \\ p + 4 (V%) ( Vl - v 2 ) \\ p 

^ cIlVViHooIKi -^|| P + c||W 2 ||co|h -V 2 ||p, 

which yields (pT)) by flB. □ 

Lemma 3.4. (i) For k G No and p > d, there exist constants c = c(k,p,d) > and 
a k G N such that for all u G W k+2 ' p (R. d ; R d ) and all X G X fc+2,p preserving the volume, 

\\V(u o X)|| fe+ljP < c||V«|| fe+1 , p (l + || VX||^ +2 + || V 2 X|p. (22) 

(ii) For p > d, there exists a constant c = c(p, d) > such that for all u G W 2,p (M. d ; M. d ) 
and X, X G X 2 '? with X - X G L p (M d ; R d ), 

\\uoX-u oX\\ p ^ c||Vu||i, P ■ \\X - X\\ p . (23) 

Proof, (i) Since X preserves the volume, we have 

||uoX|| p = \\u\\ p . (24) 

Observe that for m ^ 2 

V k (u oX) = (V m u) o X ■ (VX) m + • ■ ■ + (Vu) o X • V m X (25) 

d22]) follows by (fT8|) and (gg). 
(ii) It follows from 

uoX -uoX = ! (Vu) o (sX + (1 - s)X)) ■ (X - X)ds (26) 
Jo 

and (Jig). □ 

Lemma 3.5. For 6 No, £/ > and T := jj, there exist constants C\ = C\(p, d) > and 
c 2 = c 2 (k,p,d) > such that for any divergence free vector field u G C([— T, 0]; W k+2,p ) 
satisfying sup tg r_ T0 ] || Vwt||fe+i iP ^ U , the solution X ttS to flfy) belongs to X fc+2,p a.s., and 
for all t G [-T, 0] ' 

||VX ti0 ||oo <ci, ||V 2 X i0 |U, p < C2. (27) 



Proof. Noting that 



we have 



VX t>s = I + / {Vu r ) o Xt, r ■ VX t , r dr, 



(2f 



IIVX 



t,s I ■DC 



By Gronwall's inequality, we obtain by (1181) 



sup HVXi.Joo < exp 



||Vu r ||oodr 



T 



<; e cUT = e c =: c x . 



(29) 



On the other hand, from ( 1281) . an elementary calculation shows that 
det(VX( jS ) = exp J (V ■ u r ) o X t>r dr 
So, x i — ^ preserves the volume. Now, 

V 2 X M = J [(VV) o X t , r • (VX tjr ) 2 + (Vu r ) o X t , r • V 2 X 4 , r ]dr. 
Hence, by ((29]) and ([18]) 

l|V 2 X M || p < 0^ (llVVllpllVX^II^+llV^IUHV^II^dr 

< cUT + cU J \\V 2 X tjr \\ p dr. 

By Gronwall's inequality again we get 

||V 2 X M || P < cUTe cUT = ce c =:c 2 . 
Higher derivatives can be estimated similarly step by step. 

Lemma 3.6. For p > d and T > 0, let u,u e C([-T, 0]; W 2 > p (R d ; R d )), and X,X 
solve SDE (Tj|) with drifts u and u respectively. Then for some c = c(p, d) > and any 

te[-T,o], 

i r° 

\\X t! o -Xt,o\\p ^ ex P c SU P ll Vu *lkp 



□ 



te[-T,o] 



Proof. We have 



X t Jx) - X t Jx) 



| 1 1 ^)Ci?^« 



(u r (X t , r (x)) - u r (X t , r (x)))dr 
(u r (X ttr (x)) - u r (X t)r (x)))dr 

-s 

+ / (u r (X t ,r(x)) - u r (X t>r (x)))dr. 



(30) 



For R > 0, let Br := {x G M. d : |x| ^ i?}. By virtue of x i— > X tjt .(a;) preserving the volume 
and formula (1261). we have 



\\X t)S ~ X^ s \\ lp(Br) < / ||u r o X iir - w r o X tir \\ L p( BR) dr + / ||w r - u r \\pdr 



< SUp ||Vw r ||oo 
re[-T,0] 
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||X tjr - X ir .|| LP(i j i j)dr 



+ 



\u r — u r \\ p dr, 



(31) 



By Gronwall's inequality and (fT8l) . we get 



\X ttS — X^ s \\ LP ( Br ) ^ exp 



c sup ||Vw t ||i )P 

■ *6[-T,0] 



U-r 



u r \\ p dr. 



□ 



Letting R go to infinity gives (1311 . 

We are now in a position to prove the following local existence result. 

Theorem 3.7. For v ^ ; k £ No and p > d, there exists a constant cq = co(k,p,d) > 
independent of v such that for any uq £ jy fc+2 ' p (IR d ; M. d ) divergence free and T := 
(c || Vwollfc+i.p) -1 ; iaere zs a unique pair (u,X) with u £ C([— T, 0]; W /fc+2 ' p ) satisfying 



(32) 



X tjS (x) =x + J u r (X t;r (x))dr + V2z/ / cr r dS r , t ^ s ^ 0, 

w t = PE[(V t X ti o)(u oX tj0 )], t^O. 
Moreover, for any t £ [— T, 0] 

||VMi||fc + i iP ^ c || Vu |U+i,p- 
Proof. Set t4(a;) := «o(^)- Consider the following Picard's iteration sequence 

f X£ s (x) = x + J <(X" r (x))dr + «r(r)dB r , * ^ s < 0, 

1 ^^PEKV^KoX^)], f<0. 
Noting that 

we have by (OED and (1221) 

||V< +1 || fc+1 , p ^ E||VW(n oX™ ,X^ )|| fc+1 , p 



(33) 



(34) 



<: cE 
< c 3 E 



IVX 



t,0l|oo 



lv 2 x"|u, p )-||v( Mo oX™ 



1+ ||VX, n l|fe+3 



L t,0lloo 



k+l,p 

Vwo|| 



where £ N only depends on k and c 3 = c 3 (k,p, d) ^ 1. 
Set 

Co := C3 (l + c ^ 3 + cA)>l, 

where ci and C2 are from Lemma 13.51 Choosing U := co|| Vwo||fc+i,p and T := l/{7 in 
Lemma 13.51 we have by induction and Lemma 13.51 



sup ||Vw™|U+i, P < U, Wn £ N. 

te[-T,o] 

On the other hand, we also have by (1271) 

|K +1 || P < cEW^X^uooX^W, 



(35) 



< cE 



n||oo||no ° X^qI 



^ c||mo|| p , 

which together with ( 1351) gives the following uniform estimate: 

SUp SUp ||u"|| fc+2l p < +00. 
n6Nte[-T,0] 
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(36) 



Now by flU and pi flU), we have 



a 



n+1 



U 



m+1 1 



< cE 



SC c 



• IKo-Xto-^oO-Xfcollp 



<-<%dr, 



where c = c(p, d, [/) is independent of n, m. From this we derive that 

0. 



limsup sup || ^ 

n,m^oo te[-T,0] 



u?-u m 



t UP 



By (I36p and interpolation inequality, we further have 



lim sup sup 

n,m-*oo te[-T,0] 



\ U t ~ u T\\k+l,p 



0. 



Therefore, there is a u £ C([-T, 0]; W k+1 ' p {R d ; R d )) such that 



Taking limits for 
(1351). 



limsup sup — iZt||jfe+i )P = 0. 

n^oo t e [_T,0] 



one finds that u is a solution of 



Estimate (1531 follows from 

□ 

Remark 3.8. T7ie constant Co in ( Tff^) usually strictly greater than 1. // Co equals 1, 
then we can invoke the standard bootstrap method to obtain the global existence. This will 
be studied in the next section when the periodic boundary is considered and the viscosity 
is large enough. 

Since the existence time interval in Theorem 13.71 is independent of the viscosity u, we 
also obtain the local existence of solutions to Euler equation (j2J). Moreover, as v — > 0, the 
solution of Navier-Stokes equation converges to the solution of Euler equation as given 
below. 



Proposition 3.9. Keep the same assumptions as in Theorem 3.1. For v ^ and u £ 

W k+2,p n L 2 , let (u u ,X u ) be the solution of UM) corresponding to viscosity v and initial 
value uq. Then for any j = 0, ■ • • , k + 1, there exists c = c(k,j,p, d, ||uo||fc+2,pj ll^olb) > 
such that for all v ^ and t £ [— T, 0] 

I „ nil . . ,. ( k + 2-j In ,,1 . k+2 In 

IK - u t\\ci < C (H*I) d — p )/{ i + ~~v>. 

Proof. Note that u £ W k+2 ' p n L 2 guarantees < £ W k+2 > p R L 2 . By 

- «?) + y< + p[« • v)< - K° ■ v)«a = o, 

we have 

-c> t ||<-n t ||l = ^(A<,<-^) 2 + («-VK-(^-VK,<-« t ) 2 
= i/(A<, < - u°) 2 + ((« - «?) • V)<, < - w°) 2 



i.e., 



«c i/||A<|| 2 - ||<-u°|| 2 + ||V<| 
-<9 t |K - u°|| 2 < ^||Am^|| 2 + || Vw^Hco 



u 



1 1 2 
t \\2i 



By Gronwall's inequality and ( 1331) we obtain 

f° r /"O 

|K — m°|| 2 ^ v j || Aw^|| 2 ds • exp / ||VM^||oods 
Jt Jt 



^ cv\t\. 
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The desired estimate now follows by the Sobolev embedding(cf. [10]): for u G W k+2,p fl L 2 

IIV^mHoo 5% C/-j ) p )( i|| u\\ k+2,p\\ U \\2 > 

where a = (J + + □ 

Remark 3.10. We cannot prove a convergence rate 0{\fui) as in [12] starting from (CHJ) 
because x \— > (X^ Q (x) — X® Q (x)) does not belong to any L p -spaces. 

4. Existence of Global Solutions 

4.1. Global Existence in Two Dimensions. First of all, we recall the following Beale- 
Kato-Majda's estimate about SIOs, which can be proved as in [19, p. 117, Proposition 3.8], 
we omit the details. 

Lemma 4.1. For p > d, let u G W 2,p (M. d ; R d ) be a divergence free vector field and 
uj := cur 1m. Then, for some c = c(p, d) 

||Vu||oo < c(l + log + ||a>||i iP )(l + II^Hoo), (37) 

where log + x := maxjlogx, 0} for x > 0. 

In two dimensional case, taking the curl for the second equation in fl32l) . one finds that 

u t : = curlut := d\U 2 — d^,u t = E[tOo o X t ,o]. (38) 

From this, we clearly have 

\\u t \\p < H^olU l^p^oo. (39) 

Basing (l37j) and representation ( 1381) . we may prove the following global existence for 2D 
Navier-Stokes and Euler equations. 

Theorem 4.2. In two dimensions, for v ^ 0, k G N 0; p > 2 and u G W k+2 ' P (R 2 ; R 2 ) 
divergence free, there exists a unique global solution (u,X) to equation (3B) . 

Proof. We only need to prove the following a priori estimate: for all t G R_ 

IKHfc+2, P < e(||uo|U+2, P , k,p,t) < +oo, 
where c(||-u ||a:+2,p, k,p, t) continuously depends on its parameters. 



Following the proof of Lemma 13.51 we have 

|VX t)0 ||oo < exp 



V'Urlloodr 



(40) 



Noting that 

Vu t = E{Vu oX tfi -VX tfi ), 

we have 

\\Vco t \\ p ^ ||Vwo||p-E||VJT t ,o||oo 

and by (J3HD and flUD 

mi|i, P < lko||i,p • ^1 + exp jf HV^IIoodr 
Hence, by ([37]) dSSJ) and (jilj) 

iivtttiioo < c(i + iog + n^iii iP )(i + mi 

^ c + c / || V-Urlloodr, 

*12 



(41) 



where c = c( ||a^o|| p) - By Gronwall's inequality we obtain 

llVutHoo s= ce c| *L 
Substituting this into (j4"Ul) and (PHI gives 

II VA t( o||oo ^ e 1 1 , 
and by Calderon-Zygmund's inequality about SIOs (cf. [22J) 

||V«t||i,p < \\ut\\i, P < ||^o||i, P • (l + e c| * |eC| '' 
Moreover, 

||u t || p < cE(||VX ii0 ||oo • ||«o ° X t fl\\p) < c||«o||p • e c| * |eC|t| . 



Thus, 



IKIkp < c(||w || 2 ,p,p,t) < +00. 



Starting from fl38|) and as in Lemma I33| higher derivatives can be estimated similarly. □ 

4.2. Global Existence for Large Viscosity. In this section, we study the existence 
of global solutions for large viscosity and work on the d-dimensional torus T d = IR d /Z d . 
Let W k,p (T d , M. d ) be the R -valued Sobolev spaces on T d with vanishing mean. Instead of 
f fT2|) . we consider 



X ttS (x) = x + J u r (X tir (x))dr + V2^(B S - B t ), 



(42) 



; R d ), i.e., for to G tt, B.(u) = oj{-). 



where B is the standard Wiener process on Vt := C( 

We recall the following Bismut's formula (cf. [21 [8]). For the reader's convenience, a 
short proof is provided here. 



Theorem 4.3. For any t < and f G C^TT; R), it fcotas t/icrf 



(VE/(X t , ))(z) 



tV2z7 



=E 



f(X t>0 (x))J ( S (Vw s )oX t)S (x)-l)dS s 



In particular, for any p > d and some c = c(p, d) 



\VEf(X tfi )\\ p ^ 



'v\t\ 



\t\ ■ sup II Vm s ||i )P + 1 

se[t,o] 



(43) 



(44) 



Proof. Fix t < and 7/ G M d below and define 

1 



tV2u 



(t - s)y + / [(Vu,.) o A- t , r (a;)] • (n/)dr 



, sG [t,0]. 



Consider the Malliavin derivative of X tyS with respect to the sample path along the direc- 



tion h, i.e., 



iAAj s (x, u;) = hm — ! 1 , w G Si. 

e^O £ 



From fl42l) one sees that 



D h X t Jx) 



J [(V«r) o AV(x)] ■ D h X tjr (a;)dr + v^O) 



it - s)y 



+ / [{Vu r )oX t>r (x)) 



D h X t Jx) + 



ry_ 
t J 



dr. 
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On the other hand, we have 

VX M (x) -y = y + J [(Vu r ) o X^x)] ■ VX^(x) • ydr. 
By the uniqueness of solutions, we get 



In particular, 
Now 



VX ttS (x)-y = D h X t!S (x) + ^. 



VX t)0 {x)-y = D h X t>0 {x). 

VEf(X tj0 )-y = E[[(V/)oI t , ].VI t ,o' ! /] 
= E[(V/)oI (i0 .D,I f|0 ] 

= E[D h (foX t>0 )] 

2 (foX t , ) [ h(s)dB. 



where the last step is due to the integration by parts formula in the Malliavin calculus 
(cf. [20J) . Formula flU now follows. 

For estimation f j44l) . by Holder's inequality and Ito's isometry, the square of the right 
hand side of ( l4"3l) is controlled by 



1 



2vt 2 



E\f(X tfi (x))\ 2 E 



< ^E|/(X, (x))| 2 



J \s{Vu s )oX t)S {x)-l\ 2 ds 

t\ 3 SUp || Vttsll^, + \t 



se[t,o] 



Hence, by fllg|) 



||VE/(X t , ) 



'v\t\ 
c 



|t| • SUp || VUglloo + 1 

se[t,o] 

\t\ ■ sup ||Vu s ||i iP + 1 

se[t,o] 



□ 



The proof is complete. 

We now prove the following global existence result (see also [T6l [14] ) . 

Theorem 4.4. Let k G No and p > d, uq G W^CF*; R d ) &e divergence free and mean 
zero. Let (u,X) be the local solution of (SB) in Theorem 3.7 Then, there exist Tq = 
T (k,p,d, || Vwollfc+i.p) < and 5 = 5(k,p) > such that if v ^ <5|| Vuo||fc+i, p , then 

llVtt^Ufc+ij, ^ ||Vtt ||jfc + i s p, 

and there is a global solution to equation KSJfy . 

Proof. Let (u, X) be the local solution of ( 1321) on [— T, 0] in Theorem 13.71 where T = 
(c || Vwo||fc+i,p) _1 - Recalling the estimations in Lemma [331 and Theorem 13.71 we have for 
all t e [-T, 0] 



,o||oo <d, ||V 2 X i0 ||fc, P ^ c 2 
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(45) 



and 

||Vut||fc + i,p ^ c ||Vw |U + i,p. (46) 

Write 

u t = PE[(V%, - I)(«o o X ti0 )] + PE( Mo o X tfi ). (47) 

We separately deal with the first term and the second term. 

For the first term in (1471) . using (1451) (I46p and as in Lemma one may prove that for 
some c = c(k,p, d) and all t G [— T, 0] 

HV^o - I\\ k +2,p < c||Vmo|U + i, p • \t\. 
Using this estimate as well as (|22|) f|45l) and f|46l) . one finds that 

||VPE[(V t X t>0 - I)(«o o ^,o)]IU+i, P ^ cUVitolll+ij, • |*|- (48) 
For the second term in (1471) . by (1441) . f|46|) and Poincare's inequality, we have 



||VPE( Mo oX, i0 )|| p < -^=\\u \\ P ^ -^=:||V«o|| p . (49) 

y/u\t\ y/u\t\ 

Note that 

V 2 E(u o X ti0 ) = VE((V« ) o X f>0 ) + VE[((Vu ) o X tfi )(VX tfi - I)]. 
As above, we have 

||VE((V«o)oX ti0 )||p<^=||V« Q || p 

and 

||VE[((V«o) o I t ,o)(VI t , - I)]|| fciP < c||V Uo || 2 fc+ljP • |*|. 

So, 



||V 2 PE(u oX f , )|| P ^ c||V 2 E(n oXi, 



Continuing the above calculations we get 



\\Vu \\ p + c\\Vu \\l +hp -\t\. (50) 



||V fe + 2 PE(u o X, )|| p < -=|| V k+1 u \\ P + c\\Vu \\i +1 , p ■ \t\. (51) 

Combining (J49]) (J3DD and (|5T|) . we find 

||VPEKoX 4 , )|| fc+liP ^ - 7 £= f ||Vizo|| fclP + c||V«o||fc + i lP - |*|. (52) 

Summarizing (jUJ f|4Hl) and fl52|) yields 



llVtttUfc+ij, ^ 



+ c 4 ||Vuo||fc+i,p ■ |*| 



HVuoll 



k+l,pi 



where C3 = cs(k,p,d) and C4 = c^(k,p,d) > cq. Now, taking T = ~ 2 C4 ||vJolU+i an< ^ 
5 = 8C3C4, we have for v ^ 5||Vwo||fe+i,p 

||Vmt IU+i,p ^ ||Vw |U + i,p. 

The proof is thus finished. □ 
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5. A Large Deviation Estimate for Stochastic Particle Paths 

Let G k denote the fc-order diffeomorphism group on IR d , which is endowed with the 
locally uniform convergence topology together with its inverse for all derivatives up to k. 
Then G fc is a Polish space. Let Gq be the subspace of G fc in which each transformation 
preserves the Lebesgue measure, equivalently, 

Gq := {X G G k : det(VX) = 1}. 

Then Gq is a closed subspace of G k , therefore, a Polish space. 
It is clear that 

t » A7(-) G G fc 

is continuous by the theory of stochastic flow (cf. [IZ]). We now state a large deviation 
principle of Freidlin-Wentzell's type, which follows from the results in [U, [21] by using 
Proposition 13 .91 

Theorem 5.1. Keep all the things as in Proposition \3.9\ . For any Borel set E C 
G([-T,0];G§), we have 

- inf ICY) < liminf u\og¥(X u G E) < lim sup v log F(X U efi)< - inf /(F), 
^e£° f-+0 „_> Q y g £ 

where E° and E denotes the interior and the closure respectively in C([— T, 0]; Gq), and 
/(y) zs £/ie rate function defined by 

I(Y):=l inf r ||/i.||? 2 da, F G C([-T, 0]; G§), 

where S(h) = Y solves the following ODE: 

Y s (x)=x+ [ M°(F r (x))dr+ ! (<T r ,hr)pdr, sg[-T,0]. 

J-T J-T 

Remark 5.2. In two dimensions, the T in the above theorem can be arbitrarily large by 
Theorem \4-2\ 
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